HD-TVP-98-11 



Ferromagnetism in single-band Hubbard models with a partially 

flat band 

Andreas Mielke* 
Institut fur Theoretische Physik, Ruprecht-Karls-Universitat, 
Philosophenweg 19, D-69120 Heidelberg, F.R. Germany 

March 10, 1999 



Accepted for publication in Physical Review Letters 
Abstract 

A Hubbard model with a single, partially flat band has ferromagnetic ground states. It is shown 
that local stability of ferromagnetism implies its global stability in such a model: The model has 
only ferromagnetic ground states if there are no single spin-flip ground states. Since a single-band 
Hubbard model away from half filling describes a metal, this result may open a route to metallic 
ferromagnetism in single band Hubbard models. 

PACS-numbers: 75.10.Lp, 71. lO.Fd 

Ferromagnetism of itinerant electrons is an old and still unsolved problem in theoretical physics. One 
of the motivations to introduce the Hubbard model §\§ more than 30 years ago has been to understand 
this problem [Q]. Conceptually, the Hubbard model is very simple. It describes electrons on a lattice 
interacting via a repulsive, purely local Coulomb interaction. Originally, only a single energy band was 
considered. But it is relatively easy to generalize the model including more energy bands or further 
interactions. Despite the simplicity of the model, there are only few rigorous results concerning the 
existence of ferromagnetism. The first rigorous result is the so called Nagaoka theorem [S]. It states that 
on certain lattices, for an infinitely large Coulomb repulsion, and close to half filling, the ground state is 
ferromagnetic. But the problem is that this ferromagnetic state disappears if the repulsion has realistic 
values or if the number of electrons is varied [Q, gj, at least for usual cubic lattices in d dimensions. [||] 

Ferromagnetic ground states also occur, if one of several bands of the Hubbard model is dispersionless 
(Lieb's ferrimagnetism ^ and the flat-band ferromagnetism [Q]). Tasaki [Q] was able to show the local 
stability of ferromagnetic ground states in related models with a nearly flat band. But these models (as 
well as the flat-band ferromagnetism) are already extensions of the single band Hubbard model since they 
contain more than one band. Furthermore, these models show characteristic features of an insulator: The 
(nearly) flat band is half filled and there is a gap in the single-particle spectrum. Other extensions of the 
Hubbard model have been proposed as well, for instance an additional ferromagnetic interaction between 



the electrons [JlOp or degenerate bands with Hund's coupling between the bands []11[]. 

The question, if and under which conditions a simple, single band Hubbard model can show ferro- 
magnetism, is still open. Ferromagnetism is a strong coupling phenomenon. To obtain ferromagnetism 
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within usual mean field theory the dimensionless parameter ppU has to be large, where is the density 
of states at the Fermi level. But it should be noted that a large value of pfU is not sufficient for the exis- 
tence of ferromagnetism (counterexamples with p F = 00 and U = °° can be found in fll2|]). On the other 
hand it seems to be clear that a strong asymmetry of the band, together with a large density of states near 
the Fermi energy is a condition that favours ferromagnetism. This view is supported by the variational 
calculations by Hanisch et al. ^ and as well by the dynamical mean field analysis by Wahle et al. JlQ[] . 
In this letter I discuss the extreme limit of this situation, namely a Hubbard model with a partially flat 
band. The main results is that local stability of ferromagnetism implies its global stability in such models. 

At a first glance a single, partially flat band seems to be similar to the flat-band ferromagnetism 
mentioned above. But there are important differences: As already mentioned, the models with a nearly 
flat band have a gap in the single particle spectrum and the band is half filled This is a typical 
situation for an insulator. In the flat band case, it was possible to go away from half filling [||], but if an 
entire band is flat, single particle eigenstates may be localized. In fact, the existence of a localized basis 
was an essential point in the proofs [§]. Localized states are as well typical for an insulator. In the present 
case, the energy band is not half filled and there is no gap in the single particle spectrum. A localized 
eigen-basis of single particle states does not exist, instead the translational invariance will be important. 
These are typical properties of a metal. 

The Hamiltonian of a single band Hubbard model on a ^-dimensional translationally invariant lattice 
with periodic boundary conditions is 




c\ c (etc) creates (annihilates) an electron with spin a in a single particle eigenstate given by \\fk(x) = 

— 1/2 

N s exp(ifoc), and N s is the number of sites of the lattice, which is equal to the number of single particle 
states in the band. The interaction is usually written in the form ^Lv c v| c x| c ^l c ^T- ^ * s repulsive, U > 0. 
The Hamiltonian has the usual SU (2) spin symmetry. Furthermore, due to translational invariance, the 
total momentum is a good quantum number. The wave vectors k are elements of the first Brillouin zone 
(BZ). More precisely, k is a representative of a class of equivalent wave vectors each belonging to a 
different BZ. A statement like k + 8 = k! means that k + 8 and k' are both representatives of the same 
class of wave vectors. I assume that the single particle band is partially flat, a finite fraction of the band 
energies is degenerate. I assume that the set of degenerate band energies is situated at the bottom of 
the band, but some comments on other situations are given as well. Shifting the energy scale one can 
always take > 0. As a consequence, H > 0. Let L be the subset of wave vectors k with 8jt = 0. Nd is 
the degeneracy, i.e. the number of elements in L. In such a situation the Hamiltonian has ferromagnetic 
ground states if N e < Nd- Let 

¥of = I1 c 1tI ) ( 2 ) 

keL 

be the (only) ground state of H with S = S z = Nj/2, N e = Nd- For N e < Nd one can construct 
ferromagnetic ground states with S = S z = N e /2 by replacing the product in (Q) by a product over an 
arbitrary subset of L with N e elements. Using the SU (2) invariance of the Hamiltonian one can construct 
further ground states with S z < S = N e /2. The question is: Are there ground states with 5 < N e /2. 

\|/of has zero energy. Any other ground state with N e < Nd must have zero energy as well. Since 
both parts of the Hamiltonian, the kinetic energy and the interaction, are non-negative, the ground states 
of H are simultaneously ground states of the kinetic energy and of the interaction. This simplifies the 
situation a lot. Let us assume that H has a ground state with N e = Nd — n + m electrons (n > m) and a spin 
S = S z = N e /2 — m. Such a state can be written in the form 

y = S n _i m (a)yoF (3) 
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where 



S n l m {a)= £ a, / a a 11 V J! ^ (4) 

l\...\ m \k\...k„ j j 

OC/j ...i m -fr x ...k n are antisymmetric in the first m and in the last n indices. \|/ should be a state with spin S = S z = 
N e /2 — m. This is the case if S+\\f = where S + = Y*k c \-\ c ki, which yields Y*k a k,h...l m ;k,kz...k n = 0- Since 
the Hamiltonian is translationally invariant, the eigenstates of H are also eigenstates of the momentum 
operator. Let \|/ be a state with momentum p. This means that a/ ( ...i m ;k\ ...k„ vanishes if Y!f=\ h ~ T!j=i kj ^ 
p. Since \|/ is a ground state of the kinetic energy, CLi 1 .„i m - > k 1 ...k n nas to vanish if some indices lj are not in 
L . Furthermore, I let a/, ...i m ;h...k„ = if some indices kj are not in L . \|/ is a ground state of the interaction 
if and only if LztC/fc+ST c -&lV = for all 8. This yields a condition for a, namely 

E (-l) P ai 1 ...l m - l -k P(n + 1) +8;kp(i)---k Pi n) = ° (5) 
PeS n+ i 

for all kj, lj, 8. S n+ \ is the group of all permutations P of n + 1 objects and (— l) p denotes the sign of the 
permutation P; (— l) p = 1 if P is even and — 1 if P is odd. I let 8 = L"=i &/ — L7=/ h + Z 7 ' smce otherwise 
(Q) is trivial. Using the fact that a is antisymmetric in the last n indices, one can rewrite (||) in the form 

n+l 

£(-l)" ( ' 1) az 1 ...; m _ 1) p + £"+i^.- E ™-// / _fe j ;fc i+1 ...A n+1 ,fc 1 ...fe j _ 1 =0- (6) 

;'=1 

The sum runs now over all cyclic permutations of the indices kj. I define 

&l 2 ...l m ;k 2 ...k n = Y J a h+p.h-l m ;k l ...k n (7) 

h 

where p ^ is chosen such that 6c is not identically zero. This is possible since for some k G L, k + p 



is also in L and since a is not identically zero [13]. I put l\ = k\ +p in (S) and sum over jfci. Using the 



definition of a and the antisymmetry of a in the last n indices, one obtains 

L a *l h-lm-\ ,P+Z%\ kj-ZfZz 1 >i-h 
k\ 

H+1 

- T,(- l ) in ^ ){ ^ 2)& i 2 .j m ^,p-p+rz^ =0 - (8) 

The first term in this equation vanishes. The reason is that for each k\ there is a term in the sum over k\ 
with &i replaced by p + L"=2 kj ~ ^"j=2 h ~ ^l- Due to the antisymmetry of a in the first m indices, these 
two terms are equal up to a different sign and annihilate each other. This yields 

n+l 

£(-!)(" 2) \.j m _ up -p+jn+lk i-krJci+l-kn+i,k2-ki-i =0 W 

1=2 

which is the same condition as @ for a instead of a. Consequently \Jr = S"_ ,m 1 (ol)\^of is also a 
ground state of H. This shows that if H has a ground state with N e = Nd — n + m electrons (n > m) 
and a spin S = S z = N e /2 — m, it has also a ground state with the same number of electrons and a spin 
S = S z = N e /2 — m+l. One can now iterate this procedure to obtain finally a single spin flip state with a 
spin S = S Z = N e /2 — 1: 

Theorem - In a single band Hubbard model with a A^-fold degenerate single particle ground state 
and N e < Nd electrons local stability of ferromagnetism implies global stability. The model has only 
ferromagnetic ground states with a spin S = N e /2, if there are no ground states with a single spin flip, i.e. 
with a spin S = N e /2 — 1. 
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Remarks - 1. The existence of single spin flip ground states is relatively easy if N e = Nd- A general 
single spin flip state with momentum p is given by 

¥ = £ Uk+ P ,kcl +pl c k] \\r 0F . (1 0) 

k 

This state has a spin S = N e /2 — \ if p^O. The condition (Q) shows that this state is a ground state if and 
only if CLk+ P ,k = CLk'+pM' f° r a ^ k,k' E L. Therefore, it is possible to construct a single spin flip ground 
state with momentum p if and only if z k+p = for all k with = 0. In that situation, the single particle 



density matrix in [ ]14[ ] is reducible. Thus, for N e = Nd our result is as well a consequence of the result 
in JT4|], but due to translational invariance, the condition for the occurrence of ferromagnetism is much 
simpler. For N e < Nd it is not possible to obtain a similar (simple) condition for the existence of a single 
spin flip ground state. 

2. Let us consider situation where the degenerate single particle energy lies at the upper band edge 
of the single band, and let N e > 2N S — Nd- In this case the model has ferromagnetic ground states with 
a spin S = (2N S —N e )/2. Performing a particle-hole transformation one obtains a model that fulfils the 
conditions of the above theorem. Thus local stability of ferromagnetism implies global stability in this 
case as well. 

3. In the above derivation one uses the fact, that the Hamiltonian is translationally invariant. This is a 



natural assumption. But it is also possible to investigate a more general case. As for N e = Nd in [14], the 



proof is much more complicate and less intuitive [16] 



4. The result is true for any U > 0. U may be arbitrarily small. Therefore one may wonder whether 
this model corresponds to a strong coupling situation. This is indeed the case. The relevant dimensionless 
parameter is pf U. It is infinite in our model for any U > 0, since the Fermi level lies in the region where 
the band is flat. A partially flat band is certainly an unrealistic situation. Typically, this assumption has 
the consequence that the hopping matrix elements t^, = N^~ Y,k £ k^p(ik(x — y)) have a longer range than 
usual. On the other hand, one may hope that as in the flat band case |§] our result extends to a nearly flat 
case as long as U is not too small. This would be a realistic situation in transition metals. 

5. In most cases it is much simpler to study the stability of the ferromagnetic state with respect 
to single spin flips than the global stability. In the case of a Hubbard model with a nearly flat band, 
Tasaki [^] was able to show the stability with respect to single spin flips. In most variational treatments 
single spin flip states are used as well The variational studies of the stability of the Nagaoka state 
[|], ||] are complicated and the general single spin flip problem is too difficult to be solved completely on 
usual lattices in 1 < d < oo dimensions. In our case the situation is simpler since we already know that 
there are ferromagnetic ground states. The aim is only to show under which conditions there are other 
ground states. If there are further ground states, one should expect that small perturbations are sufficient 
to destroy ferromagnetism. This would be an instable situation. If there are no other ground states, it 
is possible that ferromagnetism is stable with respect to small perturbations of the Hamiltonian. But it 
is difficult to investigate this problem since in the present model there is no gap in the single particle 
spectrum. This is the main technical difference to the models discussed in fl. 



6. As in []14p one can generalize the above result to a situation, where the degenerate single particle 
states are not at the bottom of the band. I assume again that the A^/-fold degenerate single particle state 
has energy and that L is the subset of k with = 0. Let L < be the subset of k with < and let N < 
be the number of elements of L<. For U = and N e < 2N< +Nd the ground states are highly degenerate: 
Each single particle state with k G L< contains two electrons and the remaining N e — 2N < electrons can 
be distributed arbitrarily among the states with zero energy. If U is small one can apply degenerate 
perturbation theory. This means that among these degenerate states one has to find those with a minimal 
interaction energy. Since the system is translationally invariant the contribution from the single particle 
states with k € L < is the same for all the degenerate multi particle states. It is therefore sufficient to 
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minimize the interaction energy of the electrons in single particle states with k £ L. This is equivalent 
to the above situation, where = was the bottom of the band. If the degeneracy is lifted at first oder 
in U, the ground state is ferromagnetic with a spin S = (N e — 2_/V<)/2. Depending on N e , the spin varies 
between and Nd/2 and may be extensive if Nd is extensive. If the degeneracy is not lifted at first order, 
another ground state with a smaller spin is usually favoured at higher order in U. This argument explains 
e.g. the results presented in [15]. But this argument is perturbative and holds only for (very) small U. Is 
it possible that this ferromagnetism disappears when U becomes larger? In the moment I am not able to 
answer this question. But if the following conjecture is true, the ferromagnetism is stable for any finite 
U. 

Conjecture - Let £"05 be the smallest eigenenergy of H in the subspace of eigenstates with a spin S. 
Suppose that for some U = Uq the Hamiltonian has a degenerate ground state with a ground state energy 
£"05! = Eqs 2 and Si ^ £2- Then I claim that £"051 < £os 2 f° r ^ > ^0 if ^1 > £2- 

It is sufficient to proof this conjecture for N e < N s , since the result for larger electron numbers can 
be obtained using a particle-hole transformation. I am not aware of any (rigorous) result for the Hubbard 
model that is in contradiction to this conjecture. The physical intuition behind it is simply that if for some 
value Uo of the interaction a degeneracy occurs, one should expect that the state with a higher spin should 
have a larger kinetic energy and a smaller interaction energy, so that for higher values of U the higher spin 
is favoured. As far as I know there is no proof for this conjecture. The conjecture is trivial if S\ = N e /2, 
since the energy of a state with maximal spin (and N e < N s ) does not depend on U, whereas any other 
eigenenergy is a monotonously increasing function of U. 

Let us note that this route to ferromagnetism in single band Hubbard models naturally leads to non- 
saturated ferromagnetic states if the degeneracy is not situated at a band edge. This is similar to Lieb's 
ferrimagnetism It is even possible that a single band Hubbard model on a bipartite lattice has a 
degeneracy somewhere in the (symmetric) band. But if this degeneracy occurs in the middle of the band 
Lieb's theorem tells us that the ground state has S = 0. This is not a contradiction to the above result. 
In such a special situation one can easily see that p exists such that for each k E L , k + p is as well in 
L . Therefore, as shown in the first remark, the degeneracy is not lifted within a first order perturbational 
treatment and a second order perturbational treatment favours the singlet state, as predicted by Lieb's 
theorem. The most simple bipartite lattice, where this situation occurs, is a bipartite lattice where f^, = t 
if x and y belong to different sublattices and t = otherwise. Such a lattice is called a complete bipartite 
graph. 
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